Bond Futures Option 


We present a pricing model for bond futures option. Assuming that the bond futures price at the 
maturity of the option is lognormal, the model adopts the Black’s analytical closed-form 
solution. The futures price is taken directly from the market instead of being calculated from the 


bond futures calculator (see https://finpricing.com/lib/FxAccumulator.html). 


Let BF(t) be a price process of a given bond futures, and 7 be a payoff maturity date. The bond 


futures option with the underlying BT is a European type derivative security whose matured 


payoff at the settlement date is given by 


max {0,B(BT(T) — X)} (1) 


where £ (+1 for call, -1 for put) is the call-put index, X is the strike. Assuming that the bond 


futures price at the maturity of the option is lognormal, the prices of bond futures options at time 


zero are given by using the Black’s formula, which are 


c = P(0,T)[BF, -N(d,) - XN(d,)] (2) 


p= P(0, T)[ XN(—d, ) — BF, 7 N(-d, | (3) 


where c and p are the call price and the put price, P(0,T) is the discount factor, and d, and d, 


are given by 


_ In(BF, /X)+0°T/2 
oVT 


_ In(BF, /X)—0°T/2 
oVT ) 


d, (4) 


d, (5) 


The volatility ois defined so that oVT is the standard deviation of the logarithm of the bond 
futures price at the maturity of the option. The above pricing formulae are based on the 
assumption that the futures price is a martingale, i.e., there is no drift term in the futures price 


dynamics. 


There are closed-form solutions available for the Greek values for options described above. 


GED adopts closed-form solutions for Deltas and Gammas, which are given by 


A=e "N(d,) for calls and A=e''(N(d,)—1) for puts, (6) 
T=e" eCue for both calls and puts. (7) 
BF, -oVT 


For Vega, Theta and Rho values, GED employs numerical differentiation. Vega is calculated as 


V=P,-P, (8) 


where P, is the option price corresponding to an upper 1% shift of futures price volatility, and 


P, is the initial option price. Theta is given by 


@=P,-P, (9) 


where P, is the option value assuming one-day shift of value time while all else remains the 


same. And rho is calculated as 


tho =P. —P, (10) 


where P, is the option value assuming 10 bp parallel shift of interest rates. 


